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Abstract 



We study the phase diagram and the excitation spectra of an array of small 
Josephson junctions at / = 1/2 and arbitrary charge frustration. We find that 
the supersolid region in the phase diagram is larger than the correspondent 
region at / = and it includes two different phases. In the chiral supersolid 
(SS) charges and vortices are arranged in a checkerboard pattern on a 2 x 2 
supercell analogously to the unfrustrated case. We find a new phase, which 
we term non-chiral supersolid (NCSS), which has no corresponding phase 
at / = 0. The excitation spectra in the supersolid regions show particle like 
dispersions which is related to the defectons. The defecton condensation leads 
to superfluidity in the presence of charge ordered background. 



Typeset using REVTgX 



I. INTRODUCTION 



Josephson Junctions Arrays (JJA) are ideal model systems to study a variety of phase 
transitions induced by thermal or quantum fluctuations These latter play a major 

role if the superconducting islands are of submicron size and drive the zero temperature 
Superconductor-Insulator (SI) phase transition. The two caracteristic energy scales in the 
system are the Josephson energy J which is associated to the tunneling of Cooper pairs 
between neighboring islands and the charging energy U which is the energy cost to add 
an extra Cooper pair on a neutral island. The electrostatic energy tends to inhibit the 
Josephson tunneling: indeed a finite U produces quantum fluctuations of the phases 0's of 
the superconducting order parameter Ae**^' one each island. \i J ^ U the system is super- 
conducting since the fluctuations of the 0's are small and the system is globally coherent. 
In the limiting case J/U —>■ oo, which will be referred as the classical case, the electro- 
static energy plays no role. In the opposite limit, J <ti U, the array is a Mott insulator 
since strong quantum fluctuations of (pi prevent the system to reach coherence (Coulomb 
blockade of Cooper pairs). The SI phase transition in JJA has been studied in great detail 
both experimentally and theoretically [^]. The effect of disorder and the presence of an 
additional glass transition has been studied in Ref. [^. At the transition the system could 
be a metal with a universal value of the conductance |§|-§[. 

Frustration in a quantum JJA can be introduced either by applying a magnetic field 
or by means of a gate voltage with respect to the ground plane. The effect of the mag- 
netic frustration has been studied extensively in the classical limit ||^. The presence of the 
magnetic field induces vortices in the system and if the frustration is a rational number 
(/ = ^/^o = p/q where $ is the magnetic flux piercing each plaquette and $o = h/2e ) 
then the ground state consists of a checkerboard configuration of vortices on a g x q cell. 
A particularly interesting case is the fully frustrated situation (/ = 1/2) where the two 
degenerate ground states consist of vortex lattice with a 2 x 2 elementary supercell. The 
current corresponding to this vortex arrangement flows either clockwise or anticlockwise in 
each plaquette. We refer to this 0's configuration as a chiral (ground) state. The effect of 
quantum fluctuations for JJA at / = 1/2 has been investigated in Ref. []TU[: besides driving 
the SI transition at T = 0, quantum fluctuations affect also the superconducting regime. 
Indeed they reduce both the superconducting transition temperature and the magnitude of 
Josephson currents (though the modulus A of the order parameter is unchanged). However 
the configuration of the phases (pi and the supercurrent flow patterns are unchanged so the 
ground state is still chiral. 

Quantum effects may be modulated by means of a gate voltage to the ground V^, which 
we denote as charge frustration. Indeed the energy difference for two charge states in each 
island with n and n + 1 extra electrons may be reduced by changing Vx. Consequently 
the effects of a finite charging energy are weakened and the superconducting region in the 
phase diagram turns out to be enlarged. In the presence of charge frustration and finite 
range Coulomb interaction, the phase diagram has a rather rich structure. For instance 
different Mott insulating phases are present. They are characterized by different crystal- 
like configuration of the extra charge on each island erii which is arranged in superlattices 
whose lattice constant depends on V^. In addition new supersolid phases may appear where 
off-diagonal (superfluid) and diagonal (charge-cristalline) long range order coexist. Since 
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the original prediction by Andreev and Lifshits there were numerous works concerning 
the determination of the location and the size of these supersolid regions in the phase 
diagram |Tl|-n. Most of the theoretical investigations are restricted to the mean field 



approximation; a step beyond this, using Quantum Monte Carlo was done in Ref. JTSf for 
the JJA Model and in Ref. for the Bose-Hubbard model. The presence of the supersolid 
region is not an artifact of the mean field approximation though its size is substantially 
reduced due to the effect of fluctuations. Investigations in one dimensional systems |]T^ by 
means of Monte Carlo simulation do not show any trace of supersolid phase. 

All these investigations concerning frustration in quantum arrays considered either charge 
or magnetic frustration. The combined effect of both frustration may lead to interesting 
effects. The most striking prediction is that for certain ratios of the magnetic to charge 
frustration the JJA will be in a Quantum Hall phase [Q. The proposal is rather suggesting 
although experiments are not yet available to support it and there is not a detailed study of 
the phase diagram which would allow to located the region where the Quantum Hall phase 
should be observed. 

Motivated by these facts we make a first step in this direction by studying the phase 
diagram and the low lying excitations of a JJA in the presence of both charge and phase 
frustration. In particular we will confine ourselves to the case of fully frustrated JJA for 
which the ground state is known in the classical limit. Although we cannot make still 
statements concerning the QHE already at this stage we find new interesting results. Two 
different supersolid phases are identified which we indicate as chiral (SS) and Non-Chiral 
(NCSS) supersolid. This NCSS is a new phase and has not a corresponding stable phase at 

/ = o. 

The whole supersolid region is enlarged as compared to the case of zero magnetic field. 

The paper is organized as follows. In the next section we introduce the Quantum Phase 
Model (QPM) commonly used to study JJA and reduce it in the standard way to the XXZ 
Heisemberg model. The phase diagram of the latter model is studied at mean field level 
in section III using the 1/5* expansion. We would like to point out that the presence of 
both magnetic and charge frustrations makes the calculation non trivial already at this 
level and that it is impossible to use other standard tools as Monte Carlo simulations or 
the self-consistent harmonic approximation. In the case of Monte Carlo there are sign 
problems while in the harmonic approximation the discreteness of charges which is essential 
to describe supersolids is not accounted for. In section IV we study the spectrum of the 
low-lying excitations in all the regions of the phase diagram. Finally, section V is devoted 
to the conclusions. 



II. THE MODEL 

The physics of JJA made of submicron size junctions is usually described in terms of 
quantum dynamics of the phases (pi (at low temperatures the fluctuations of the modulus 
A are unimportant). Since in nanofabricated samples there are no Ohmic currents between 
the islands and also quasiparticle tunneling is negligible all the relevant physics is captured 
by the QPM 

Hqp = "^{ni - 71:,) Uij {rij - n^) - J ^ cos {(pi - (pj - Aij) , (1) 

id <i,j> 
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where Ui and 0j are canonically conjugated with, [(pi, rij] = 2ei Sij. The Coulomb interaction 
is described by the matrix Uij = Ae^C^j^ where Cij is the inverse of the capacitance matrix. 
The external voltage enters via the induced charge en^ and fixes the average charge on 
each island. A perpendicular magnetic field with vector potential A enters the QPM in 
the standard way through Aij = || // A-dl and gives the magnetic frustration parameter 
/ = 2^ ^ij where the summation runs on an elementary plaquette. In the classical limit, 
J/U ^ oo, JJA's are a physical realization of the two dimensional XY model. 

We will consider very large onsite Coulomb interaction and very low temperatures where 
few charge states are important. In this limit the relevant physics is captured by considering 
only two charge states per island and the QPM is equivalent to a spin-1/2 Heisenberg model 

in 

Hs = -hJ2S^ + J2S^ U,, S]-JY: (e'-""^ S+ Sj + e^^- Sr) , (2) 

'i- hj (id) 

where the operators S^, Sj' are the usual 5^(2) operators, S'f being related to the charge 
on each island (rij = + 1), and the raising and lowering Sf operators corresponding to the 
"creation" and "annihilation" operators 6^**^^ of the QPM. The "external" field h is defined 
as h = {rix — 1/2) Uij. A particulary interesting case is when is close to a half integer 
and the two charge states are almost degenerete. The various magnetic orderings in the XXZ 
hamiltonian correspond to the different phases in the QPM: long range order in (S^) and 

{Sy) i\/{Sff + {Syf ^ 0) indicates superfluidity in the QPM, while long range order in 
{S^) (I (5*^) 1 7^ 0) describes order in the charge conflguration. 

In order to treat the / = 1/2 case we consider a square lattice and use the gauge 
A = Hyk. This gives = tt ^ sign(xj — Xi) where the coordinates of the neighbouring 
sites i and j appear and a is the lattice spacing. This is illustrated in Fig.(|l]) where the 
dashed lines indicate the antiferromagnetic bond {Aij = vr) and the continous lines the 
ferromagnetic bonds {Aij = 0). 

In the classical limit there are two degenerate 2x2 periodic ground states: in the QPM 
language they correspond to two characteristic (pi patterns (one of them is shown in Fig. |l]) 
which determine a checkerboard arrangement of vortices and antivortices. 

In the XXZ version this corresponds to (S'f ) and {Sf) reproducing the above (p pattern, 

arctan ((5'f)/(5'f )) = (pi and and \J {Sf)"^ + {Sf)"^ = 1. Charge are completely delocalized 
which corresponds to (Sf ) = 0. 

In the rest of the paper we consider the Coulomb interaction only between nearest neigh- 
bour (f/i) and next nearest neighbour (f/2) sites. 

III. PHASE DIAGRAM 

In this section we obtain the phase diagram for / = 1/2 in the presence of arbitrary charge 
frustration. Previous work already established the phase diagram of the spin hamiltonian 
in zero magnetic fleld |[TT|- [13| , p!6| . The chiral conflguration of the classical ground states 
suggests the introduction of four sublattices indicated by the index I = a, /?, 7 5 (Fig. p. 
Each site is parametrized by / and an additional index p which runs within each sublattice. 
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We now focus on the XXZ model which we study using the l/5'-expansion (S is the 
modulus of the spin vector). First the quantization axis in each sublattice is rotated to 
align the spins along the direction of the relative magnetization. We obtain the rotated 
hamiltonian Rif^R^^, where 

R = H n R^piR^pi 5 (3) 

P l=a, /3, 7 (5 

with R^pi = e'^'-^p^ ,R^pi = e'^'^pK Then we perform the 1/S expansion and determine 
the actual values of the eight angles (pi and 6i by minimizing the S ^ oo limit ifoo of 
the rotated hamiltonian. The above procedure allows to characterize ground states whose 
properties are uniform within each sublattice. The order parameters will be sinOi for global 
phase coherence and cosOi for charge ordering. The (pi configuration will determine the 
supercurrent configuration in the flux phases. 

The 1/5* expansion of the rotated hamiltonian can be carried on sistematically by using 

the Holstein-Primakoff transformation S^i = V^S ttpi ^Jl — Upi/S , S~i = (5*^) , Spi = 

Upi — S , where api , a}pi are bosons, and Upi = a^pfipi These operators describe excitations 
around the S ^ oo ground state (and not particles in the QPM). We obtain 

YiHs^^ = H^ + Hsw + O i^fs) . (4) 



where H^o is of order S*^, Hsw describes low energy fluctuations around the S ^ oo ground 
state and turns out to be of order S. The S oo ground state properties are obtained by 
minimizing ifoo which is a sum over p of identical terms of the form 



The 
the 



-f^oo = ^ XI -^COs6'i + ^C0S6'; X COS 6^^ + ^ COS 6*; X COsOm 
I m=nn{l) m,=n,nn(l) 

- ^ sin 9i e'^'"' sin 9m COS {(pi 

m=nn{l) 

; minimization of the ifoo leads to eight equations in the variables (pi and 9i defined on 
plaquette of Fig.(|l|.b). Four of these equations can be solved analitically with the result 



— ' 4 ^ — ' 

m=nn{l) m=nnn{l) 

(5) 



tan (0^ — 
tan [(ps — 



'iQaW - (62 - 1) (a2 - 1 
(62 - 1) {a? + 1) ' 
16a262 - (62 - 1) (a2 - 1) 
(a2 - 1) (62 + 1) 



where 



JlQaW - (62 - 1) (a2 - 1) 

tan (0, - = -j^^^Y)W^vrm^) • 



sin 9 a sin 9r 
sm 9 s sm t/^ 



sin 9s ' sin 9^ 

We now look for solutions with a = 6 = 1 which correpond to phases with checkerboard 
simmetry. The resulting configuration of the (pi is the same as in the classical limit, so 
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the supercurrent has the same chiral pattern, the ground state is chiral as in the classical 



situation. This has been noticed in the case of zero external charge in Ref |10[. We show 
that for non zero charge frustration the above conclusion holds true for the superfluid phase 
and for the chiral SS supersolid. 

The remaining equations for di reduce to the ones of the / = case with the substitution 
J ^ the solutions then read 



cos 9 



h 



2 (Ui + U2 + V2J 



(7) 



for the superfluid (SF) and 



cos6'q, 
cos 9fj 



4w 



I [V + y/v'^ — Aw 



(8) 



the chiral SS supersolid. We defined w = {h - hs)/ {2U2K,) X = 2V (f^i - U2) - 2J2 , 



1- [V2J/ {Ui-U2)\ (l + w). 
By compairing the energies relative to the above solutions we obtain the boundaries 
between the phases with ceckerboard simmetry 

Paramagnetic - canted state 



h = ±2 (Ui + U2 + V2J 



(9) 



Canted state - ^ or | lobe 



h = ± 



U1 + U2 + V2J + Mui + U2 + V2j) (ui + U2- 3V2J) 



(10) 



Canted - chiral supersolid state 

h = 2(Ui + U2 + V2J 
Chiral supersolid state - Neel state 



U1-U2- V2J) 
^ (f/i -U2 + V2J) ' 



(11) 



/i = 2V(f/i-t/2) -2J2 



Chiral supersolid - 7 or | lobes 



h = ± 



2U2 + hs~ V(2f^2 + hsY - hs' - 8U2 (Ui - U2 



(12) 



(13) 



The spin configurations corresponding to the insulating, superfluid and supersolid regions 
are shown in Figs. ( 0). In the Mott lobes there is no projection of the spin on the xy plane and 
the three spin configurations in Fig.( |^.a) correspond to filling 1, 1/2 and 3/4. In Fig.( ^.b) 



6 



and Fig.( Be) the superfluid and supersolid regions are represented. The fact that the 



configuration is not affected by quantum fluctuations has aheady been noticed in Ref. |T0 
for zero external charge. We extend this result to the SF and SS phases resulting at non 
zero charge frustration. However we discuss below the NCSS solution which does not show 
the above supercurrent pattern. 

The NCSS ground state is found if we look for solutions with no ceckerboard simmetry. 
A non-ceckerboard phase was also found at / = by Bruder et al [] and was named SS2 
supersolid. We flnd two degenerate NCSS conflgurations namely (pa = 4>f3 = (p-y = ^ {.4>s is 
meaningless since sinO^ turns out to vanish) and (pa = 0, (pfs = (ps = tt, (sin6'^ = 0). We 
obtained the regions of the phase diagram where the NCSS solution is stable by minimizing 
numerically the energy (^). The resulting pattern of the supercurrent is different from the 
classical chiral one. We point out that in this case the equations for the angles 6i cannot be 
obtained from the / = equations by simply rescaling J ^ J / a/2. 

The solutions are characterized by the sinO order parameter vanishing on one of the 
four sublattices. On this particular sublattice cosO = 1 so the charge is well deflned. Thus 
the NCSS ground state describes a supersolid, since charge-crystalline and superfluid order 
coexist, but the chiral supercurrent pattern is lost. The phase boundaries obtained numeri- 
cally and the conflgurations of the angles (pi and di of the NCSS ground state are shown in 
Fig.(|.d). 

All the phase boundaries obtained in this section are summarized in the phase diagram 
shown in Fig.@. The whole supersolid region (SS and NCSS phases) is enlarged compared 
to the / = case (SSI and SS2 phases). At / = the tip of the lobe (which coincides with 
the extension of the supersolid phase in the hard core limit) will correspond to J/f/i = 0.45. 
A blow-up of the NCSS region is shown in Fig.(P and compared with the rescaled / = 
phase diagram. This flgure emphasizes that although the whole supersolid region is obtained 
by rescaling the / = phase diagram, the individual phases are not (the rescaled SS2-SS1 
phase boundary is shown by crosses). 



IV. EXCITATION SPECTRA 

The excitation spectra can be obtained calculating the eigenmodes of Hsw- As pointed 
out in the previous section all the terms up to O (S) are retained in the 1/S'— expansion; this, 
in turns, corresponds to retain products at most bilinear of the creation and annihilation 
operators. If we regarded / as a colour index, Hsw describes a system of four kinds of 
interacting /-bosons deflned on a lattice whose sites are labeled by the index p. Then we 
Fourier transform with respect to p and Hsw reduces to a sum of single mode hamiltonians 

Hsw = y^,Hk , 

k 

y^k +^k +^k )^ka+\^k +^k +^k + 

[yk 0.k,aO'krf + % «fc,7«fc,a ) + [^k (^k,a(^k,f3 + % (^k,l3(^k,a ) + 
( t , t \ ^ I (7,5) t , (t,^)* t A , 

[yl 'alpak,5 + vl alsak,i3) + [vk^ '' 'al^ak,5 + vl al sak,y) + 
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(^Q^fc (^k,a(^k,j + % ak,jak,a) + [Qk (^k,a(^k,f3 + (Ik (^k,l3ak,a) + 



(^k,i3(^k,s + (Ik (^k,pak,s ) + [Qk (^kn(^k,6 + (Ik (^krf(^k,s ) + 



»,S) t t I V 



0'k,aO'k,5 J + ^^Q^fe C'k,pC'k,'y + ?fc Q'k,'yO'k,l3 



(14) 



The coefficients in the hamiltonian ([14D (reported in the Appendix A) depend on the 
angles 9i and 0/ introduced in the section III. The fc-sum is restricted to half of the Brillouin 
zone because of the doubling of the lattice constant due to the magnetic frustration. Linear 
contributions in the operators ,ak,i vanish for the {6i , 0;} ground state configuration. 

The hamiltonian Hsw is diagonalized using an algebraic technique p3| , p5| briefly de- 
scribed below: we introduce the following operators 



Xw 
X/' 



Clk,lCl-k,l' + Clk,l'0'-k,l , 

(^k,i'(^~k,i + (^k,fi-k,i' = iXii>)^ 

C'k,l'C'k,l + C'-k,l'C'-k,l ! 

C'k,lC'k,l' + C^^-k,lC'-k,l' = i^Xi'' ^ 

nk,i + n_k,i + 1 , 



(15) 



which obey the following commutation rules 



[Xw , X. 
[Xw , X 
[Xw , X. 

[x'',x, 

[X/',X 



mm' 
mm' 



[X'^X™™'] =0 

; di> + Xi Ov + Xp Oi + All di 

XX m' I -y- r m' 
ImOl' + AiimOl , 

Xlm' X I' vl'nT-' X I 
Om — A Om , 

V I' xm! vl"^' X 
Am Oi — A Om , 



(16) 



where we have omitted the /c-index in the X's and {1,1' ,m,m'} are colour indeces. The 36 
operators {X;/',X"',X/ ,Xi/',X/} form a base for the non-compact symplectic Lie algebra 
sp{8)k The diagonal operators X/ generate the Cartan sublagebra of sp{8)k whose 

dimension is 4 (equal to the rank of the algebra). The other off-diagonal operators, in 
number of 32, are the non-Cartan generators of the algebra. By using the definitions of 
Eq. (pISD , one can see that Hk belongs to sp{8)k since it can be written as linear combination 
of a subset of its generators 



E + E Ar^'^ X/' + Y: Rw'^'^ Xw 

I l^l' li^V 



(17) 



where Rw^^^ = Ri'i^''^* and Dw^''^ = Dm^^^* . (The coefficients in the expression ((T^) are 
listed in the Appendix B). Eq.(|T^ suggests that Hsw posseses A = 0A;Sp(8)^ as dynamical 
algebra, therefore we can diagonalize the Hamiltonian using the fundamental faithful ( i.e. 
preserving the commutation rules) Irreducible Representation (IRR) of A. It is worthwhile 
to notice that because of the non compacteness of sp{8), every finite dimensional IRR of such 



8 



algebra is not hermitian [^]. Despite of this fact we can use IRR of A in order to obtain the 
correct eigenvalues of Hsw (see Appendix A). Such eigenvalues still depend on 9i and 0^: the 
actual spectra in the various phases are worked out by specifying the correspondent S ^ oo 
ground state {6i 0;} configuration for each phase. The main results of this procedure are 
described below for the various regions of the phase diagram. 

Insulating phases: Since the superfluid order parameter is zero there is no effect due 
to magnetic frustration. The Neel and insulating ^ (and |) phases are charge-modulated 
solids. The lowest lying excitations are gapped and particle like. They correspond to 
particle - hole excitations. In Fig. ( ^a) and Fig. ( ^.b) we show the spectrum for the 1/2 
and 3/4 lobes. The Neel solid is charge- modulated along both x and y directions and the 
gapped branches are shown in Fig. ( ^.a). The four branches in the 3/4 lobe come from 
the more complicated structure of the elementary cell. In this state the transverse phonon 
like excitation, characteristic of the diagonal long range order along rows and/or columns 
(where the charge is uniform) In addition we find the particle like spectrum (with positive 
curvature at small k) which reflects charge modulation. This is shown in Fig. ( ^.b) where 
the the dashed branch represents the phonon like excitation. 

The Mott phases are not modulated in any direction, then the low lying excitations are 
phonon like . 

SF phase: Because of the doubling of the elementary cell of the lattice due to the 
magnetic field, a gapped branch is present in addition to a gapless mode which linearly 
depend on k at small k, see Fig. ( This was already discussed in the classical case by 
Ariosa et al. 0). Long range diagonal order existing along any lattice direction induces 
transverse phonon like excitations. 

Supersolid SS phase: In the SS phase both diagonal and off-diagonal order modulation 
are present. In this state the in-phase density fluctuations are coupled to the off-diagonal 
fluctuations and decreases the sound velocity because it reduces the superfluid density. The 
supersolid is characterized also by a modulation of the charge. Such modulation of the 
number of particles can be considered as defects of the lattice and their oscillations can be 
regarded as a fluctuation of the diagonal order on each site The localized quasi-particle 
associated to the collective oscillations of the defects were termed as defectons The 



gapped branch reveals the particle like excitations consistently with this physical picture. 
In Fig.( 1^) the acustic branch is reminiscent of the superfluid order while the gapped branch 
with positive curvature is the excitation spectrum for the defectons. 

NCSS phase: As in the previous case, the off-diagonal order is coupled with the density 
waves, but in this state there are four branches. The gapless one, reflects the superfluid 
nature of this state. Since the mean fleld superfluid order parameter vanishes in a lattice 
site of the 2x2 plaquette the sound velocity is strongly decreased. Two of the three gapped 
curves take into account the defects of the lattice, the other one means that in the system 
there are transverse phase-phonon excitations too. This is summarized in in Fig.( j^). 

An important issue which can be addressed by studying the excitation spectrum is the 
determination of the dynamical critical exponent at the various phase boundaries. It turns 
out that the fully frustration does not change the critical exponents at the phase boundaries. 
At the S'S'-Neel solid transition the gap of the lower branch of the solid vanishes as k"^ giving 



a critical exponent ^ = 2 [!(:]. At the SF- S'S' supersolid transition the critical mode is 



k = (tt, tt) and the roton minimum is at this wave vector. At the phase boundary the roton 
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gap disappears linearly in k giving z = 1. 



V. CONCLUSIONS 

In this paper we considered the properties of a fully frustrated quantum Josephson 
Junction array in the presence of arbitrary charge frustration. We determined the mean 
field phase diagram at T = as well as the low lying excitation spectrum using a l/S* 
expansion for the equivalent XXZ model. 

At / = 1/2 two kinds of supersolid phases, in addition to the ordinary superfluid and 
insulating phases, are present. Besides the chiral supersolid SS which has an analog in the 
absence of magnetic frustration we find the new NCSS supersolid phase which is non chiral. 
In the case of the superfluid and the SS supersolid the main effect of quantum fluctuation and 
charge frustration is to lower in magnitude the superfluid order parameter. The supercurrent 
pattern induced by the external magnetic field is however unaffected. In our mean field 
analysis this is reflected in a rescaling of J. In other words the magnetic frustration fixes 
the current distribution while the charge frustration is responsible for the reduction of the 
magnitude of the superfluid order parameter. In the non chiral supersolid NCSS diagonal 
and off-diagonal long range order coexist in a non trivial way and the combined effect of 
charge and magnetic frustration cannot be separated as in the other phases. The NCSS 
ground state has no checkerboard simmetry since the superfluid order parameter vanishes 
in one of the four sites of each plaquette and two flux quanta may be accomodated in four 
neighboring plaquettes. There is non chiral symmetry in this phase. The NCSS phase exists 
only in the fully frustrated case and has no analog at / = 0. The whole (SS plus NCSS) 
supersolid region in the phase diagram is enlarged compared to the unfrustrated case: the 
supersolid state better adjust to the periodicity (on a 2 x 2 plaquettes elementary cell) 
induced by magnetic frustration than the SF phase and it is more favoured than in the 
/ = case. This may be important for the experiments because the combination of charge 
and magnetic frustration may help in detecting the supersolid phase. 

We also determined the low lying excitation spectrum of the system. Due to the combined 
presence of magnetic and charge frustration the excitation spectra become more structured. 
They can be revealed, for instance, by studying the anomalies in the I — V characteristics 



when Andreev current is injected into the array is coupled from a normal metal electrode 

We are currently investigating different rational values of the magnetic frustration in 
order to study the possibility of a Quantum Hall phase as predicted in Ref. 
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APPENDIX A: 



In this Appendix we list the coefficients of the hamiltonian (|Tip 

^(/,m) _^ 2^6*^'" sin 9i sin 9^ cos zi^m — cos 9i cos 9^ , 

The coefficients of the off diagonal operators are 

+ 1) cos Zl^rn - i (cOS 9i + COS 9m) Sm Zl^rn] + 

— 1) COS zi^m + i (cos 9i + COS 6^^) sin Zi^m] + 

(A2) 

where I and n indicate n.n. and n.n.n. sites respectively. In the preceding expressions 
the lattice spacing vector a/ ^ has both normahzed components and zi^m are the superfluid 
order parameter's out of phase — (pm) on nearest neibourgh sites. In the present case the 
coefficents ( |A2|) are C-numbers and they reduce to real ones in the zero magnetic field case 
only. 



4''"^ = cos (k ■ az,^) {te^^"" [(cos^; cos^, 
^sin6'isin6'm} , 

= COS (k ■ a;,^) {te^^"" [(cos^, cos^, 
^sin6';sin6'm} , 



APPENDIX B: 

The IRR faithful representation of sp{9>) we used consists in mapping the hamiltonian 
Hsw in the matrix Jvi (Hsw)- Since that sp{S) is a non compact algebra turns tu be 
not hermitian and it has the following structure (we use the same notation as in Ref. p8[| ) 

f V TZ \ 

[-n -V j 



(Bl) 



where the tilde indicates the reflection in the minor diagonal {TZ = TZ) and T? is hermitian. 
The matrix elements of "D are 



Dn^'^ = 






+ 4-" 


+ h cos 9a^ 


D,,^'^ = 


1 (MP) 






+ h COS 9/3^ 


Dss^'^ = 








+ h cos 9/3^ 


Dj'^ = 


2 l^fc 






+ h cos 9sj 



(B2) 
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2 ^ ' 2 
The matrix TZ has the following structure 



„.(7,5) 



(B3) 



I' Ru Rl3 Rl2 ^ 

-R24 -R23 R21 

-R34 -R32 -R31 

V -R43 -R42 -R41 / 



(B4) 



and its elements are given below 



p (k) ^ l^(a,/3) D (k) ^ l„(a,7) p (k) 



R2,^'^ = ^qk 



(13,1) 



R2,^'^ = ^Ik 



tl4 



^Qk , 



(7,5) 



(B5) 



The diagonalization of Hsw is equivalent to an inner automorphism of the algebra on 
itself. In other words, we can define the unitary operator U = Hfc f^fc? with 



(B6) 



where the antihermitian operator Gk is a linear combination of the non-Cartan generators 
of the algebra 



Gk = Yl [^r,s {Xrs - + ^/./ (X/ - X/)] 



(B7) 



The rotation of the hamiltonian trough U defines an inner automorphism of the algebra 
generalizing the Bogoliubov transformation 

VHsw"^-^ = E exp {ad Gu) Hu , 



exp {ad Gk) Hk = UkHkUf 



-1 



Hk + ^—f [■■■{Gk, [Gk, Hk 



n 



(B8) 



For the closure property of the Lie algebras the equation ( P8|) still produce an element of 
A, however we can fix the V'rs's and the ipr^^s in such a way to put to zero the coefficients 
of the off-diagonal part of the hamiltonian. As a final result Hsw becomes proportional to 
the generators of the Cartan subalgebra of A 



IJHsw^ ^ = HdiaG. = Uk^^^Xr 



(B9) 



r=l 



In the present case Hsw should be diagonal after having solved 16 coupled non linear equa- 
tions. Being interested in the eigenvalues only, we are allowed to use the faithful representa- 
tion of A to find the spectrum of the hamiltonian diagonalizing the matrix A4 {Hsw) instead 
of Hsw as operator using the equation (P8|). These operations are equivalent because the 
diagonalization procedure involves commutators between the elements of the algebra only. 
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FIGURES 



FIG. 1. Due to the non trivial periodicity due to the magnetic frustration a 2 x 2 ceh should 
be considered. The labels for the sites as used in the paper are here reported a). Using periodicity 
boundary conditions, as shown in a), one can reduce the problem defined on the plaquette drawn 
in b). 

FIG. 2. The 3-components spin vectors we used are characterized by the usual Euler angles 
9 and (j). a)In the insulating phases the order parameter vanishes and only the 6 configuration 
is important. b)The mean field order parameter's phase configuration is unaffected by magnetic 
frustration, c) The 9 and (f) configurations in the supersolid SS. As in the superfluid state, the 
magnetic frustration preserves the chiral order of the ground state. d)The two degenerate ground 
states of the NCSS. The circle indicates the lattice site where the phase of the superfiuid order 
parameter is not defined. 

FIG. 3. The phase diagram for the fully frustrated JJA. The n.n.n. Coulomb interaction is 
?72 = 0.1 Ui. 

FIG. 4. The region of the phase diagram containing the NCSS is shown in detail. In addition 
we report (with crosses) the phase boundary that should have the chiral supersolid SS2 by rescaling 
J as discussed in the text. 

FIG. 5. a)The excitation spectrum of the Neel insulator, at J = Ui/A, U2 = SUi/8 and 
h = 2J7i/3. The behaviour at small k reveals the particle like nature of the excitations. In b) we 
show the four branches of the "3/4" insulating phase {1/2 = O.lUi, h = 2U\ and J = O.lUi). The 
two lower curves have been rescaled by a factor 10. 

FIG. 6. The acustic and transverse phonon like branches in the superfluid are shown for 
J = lUi, U2 = 0.1 Ui, h = 2.2 Ui, ky = 0. 

FIG. 7. In this figure we show the dispersion relations in supersolid SS (J = 0.35 C/i, 
U2 = O.lUi, h = 1.7Ui, ky = 0). The sound velocity of the acustic branch is reduced com- 
pared to superfiuid case. The gapped branche has a positive curvature at small k characteristic of 
the particle like nature of the excitations. 

FIG. 8. The excitation spectrum of the NCSS phase ( J = 0.2 Ui, U2 =0.1Ui, h = 1.848 Ui, 
ky = 0). The two lower curves have been rescaled by a factor 5. 

FIG. 9. At the phase boundary between the Neel insulator and Supersolid SS the excitation 
spectrum vanishes as k^. The transition is signaled by softening of the acustic branch proper of 
the supersolid {J = O.lUi, U2 = O.lUi, ky = 0). 
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FIG. 10. At Superfluid-Supersolid SS phase transition the roton mode vanishes as k 
(J = 0.5 Ui, U2 = 0.1 Ui, = ky). 
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